Superconducting charge-ordered states in cuprates 
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Motivated by recent neutron scattering and scanning tunneling microscopy (STM) experiments on 
cuprate superconductors, we discuss charge-ordered states, in particular with two-dimensional charge 
modulation patterns, co-existing with superconductivity. We extend previous studies of a large- N 
mean-field formulation of the t — J model. In addition to bond-centered superconducting stripe 
states at low doping, we find checkerboard-modulated superconducting states which are favorable in 
an intermediate doping interval. We also analyze the energy dependence of the Fourier component 
of the local density of states at the ordering wavevector for several possible modulation patterns, 
and compare with STM results. 
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I. INTRODUCTION 

A series of recent experiments have highlighted the 
importance of spin and charge ordering tendencies in 
the cuprate superconductors. Static stripe., order has 
been established in Nd-doped L^-zSr^CuOjil which ap- 
pears to co-exist with superconductivity at very low tem- 
peratures. Recently, superconducting samples of un- 
derdoped YBa2CusO6.350 and nearly optimally doped 
B^S^CaC^Os+iu with T c up to 90 K have been found 
to show signatures of charge ordering. In addition, a 
variety of compounds display strong dynamic spin and 
charge fluctuations, which are enhanced by the applica- 
tion of moderate magnetic fieldsoEru. This emphasizes 
that even materials, where no static order can be ob- 
served in the absence of an external field, are in close 
proximity to a critical point where spin and/or charge 
order is established in the superconducting state. The- 
ories based on this-.assumptionLll3'u successfully explain 
a number of NMREj and neutron scatteringBtrU experi- 
ments. 

In this paper, we will focus on superconducting states 
with static charge order, but dynamic spin fluctuations - 
this appears to be realized in the experiments of Rcf. |^. 
Such states are found upon doaiiig paramagnetic Mott 
insulators on the square latticeE3ii3. The undoped quan- 
tum paramagnet has broken translational symmetry as- 
sociated with spontaneous bond charge (or spin-Peierls) 
order; at small carrier concentration, 5, this order per- 
sists and co-exists with anisotropic superconductivity; a 
d-wave superconductor with full square lattice symme- 
try appears above a critical S. ["Charge order" is de- 
fined very generally as spatial modulation in any SU(2)- 
invariant observables, such as local density of states 
(LDOS) per site, or kinetic exchange energy per lattice 
bond; the modulation in the total site charge density can 
be small due to long-range Coulomb interactions.] A 
particular feature of the superconducting charge-ordered 
states found in Refs. [L2[ 1^ is that the modulation is 
bond-centered and its real-space period p always takes 
even integer values (in units of the Cu lattice spacing); 
this is in contrast to the continuous doping evolution of 



the ordering wavevector, l/p, usually assumed in the so- 
called "Yamada plot" til. Interestingly, the experimental 
results of Refs. |, ||, || appear to be remarkably well de- 
scribed by the states proposed inJleis. [12], [b^, as has 
also been discussed in recent workEj'tj which appeared 
while this paper was being completed: Neutron scat- 
tering on underdoped YBa 2 Cu3O6.350 shows charge or- 
der with a real-space period of 8 lattice sites, whereas 
Bi 2 Sr2CaCuoD8+<5 close to optimal doping displays p = 4 
modulationEla; furthermore the STM data indicate an 
even-period modulation not only in the site charge den- 
sity, but also in the bond-kinetic energy and perhaps the 
bond pairing amplitudeti 

We note that numerical studies of the t — J modeltll 
have observed bulk charge order with period p = 4 at dop- 
ing level 1/8; and paired hole states with different types 
of charge order in both insulators and superconductors 
have been discussed elsewhereE3il3o. 

In the past, most theoretical work has been focussed 
on states with one-dimensional (Id) charge modulation, 
often referred, to as stripes. However, recent STM 
experimentstra indicate LDOS modulations in both x 
and y directions in a single Cu02 plane (although a 
small anisotropy is observed). Theoretically, charge 
density wjave (CDW) fluctuations are expected in both 
directional^ on the disordered side of the charge or- 
dering transition (if the system has no intrinsic lattice 
anisotropy). Moving to the ordered side, it depends on 
microscopies whether CDW order in one or in two direc- 
tions condenses, leading to stripe- like or two-dimensional 
(2d) modulations, respectively. 

The purpose of this paper is twofold: In Sec. Q, we re- 
examine the mean- field theory of Refs. [j^, [ll| to investi- 
gate the possible existence of and the mechanism leading 
to superconducting charge-ordered states with 2d charge 
modulation. In Sec. HI, we turn to a detailed discussion 



of the charge modulation pattern, by calculating the en- 
ergy dependence of the LDOS Fourier component at the 
ordering wavevector and comparing it with the measure- 
ments of Ref. S. 
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II. MEAN-FIELD THEORY 

We start with re-analyzing the largc-A theory of 
Refs. [12], [l3], which provides a microscopic description of 
doping mabile charge carriers into a paramagnetic Mott 
insulatoinJ. We consider an extended t — J Hamiltonian 
for fermions, c,- a , on the sites, i, of a square lattice with 
spin a = 1 . . . 27V (N = 1 is the physical value): 
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is the on-site charge density, and the 



spin operators are fermion bilinears times the trace- 
less generators of Sp(27V). For most of the following, 
the fermion hopping, £y, and exchange, Jy, will be re- 
stricted to nearest-neighbor terms, t and J; for the de- 
tailed comparison with the Ref. ^ we will introduce 2nd 
neighbor hopping, t' . The electronic Coulomb interac- 
tion is represented by the on-site constraint rii < N and 
the off-site repulsive interactions Vij — V/|r^-|. The av- 
erage doping (5 is fixed by = NN„(1 — S), where 
N s is the number of lattice sites. To proceed, we rep- 
resent the spins by auxiliary fermions fi a and the holes 
by spinless bosons 6j, such that the physical electrons 
Cia = bjfia, and the necessary Hilbert space constraint is 
implemented by Lagrange multipliers Aj . Via a Hubbard- 
Stratonovich decoupling of the antifcrromagnetic interac- 
tion we introduce link fields Qij , defined on the bonds of 
the square lattice. After taking the limit N — > 00, the 



slave bosons 6, condense, (bi) = VNbi, the Qij and A, 
take static saddle-point values, and we are left with a 
bilinear Hamiltonian which can be diagonalized by a Bo- 
goliubov transformation. At the saddle point, the slave 
boson amplitudes fulfill J2i b\ = N S S, and the link fields 
are given by NQy = (J a0 fjj^) , where J afi is the anti- 
symmetric Sp(2A) tensor; for further details see Ref. |l3| . 

Various ground states obtained from the numerical so- 
lution of the above mean-field equations have been dis- 
cussed in Refs. |l^, At 6 = the ground state is a fully 
dimerized, insulating spin-Peierls state, i.e., it has period- 
2 bond charge order. At low doping 6, the bare large-iV 
t — J model tends to phase separation, and the inclusion 
of moderate Coulomb interaction, V, leads to the for- 
mation of bond-centered, superconducting stripes with a 
Id charge modulation. Large doping destroys charge or- 
der and leads to a pure <i-wave superconducting ground 
state, in this regime the large- N approach reduces to the 
usual BCS mean-field theory with renormalized hopping 
matrix elements. 

Motivated by the STM results of Refs. |, | we have 
searched for additional saddle-point solutions with 2d 
charge modulation; we have restricted the attention to 
states in which the charge distribution respects the 90 de- 
gree rotation symmetry of the lattice. Interestingly, there 



are several such saddle points which were overlooked in 
Ref. [l3|. In most of the low and intermediate doping re- 
gion the states with Id and 2d modulation are close in 
energy; at small doping the Id stripe-like states are pre- 
ferred, whereas the 2d checkerboard-likes states are lower 
in energy in a certain interval of intermediate doping and 
Coulomb repulsion. 

We have therefore concentrated on states with a real- 
spacg periodicity p = 4, leading to a unit cell of 4 x 4 
sitesd. (Recall that bond-centered stripes with the spa- 
tial period pinned to four sites were found over a rather 
large range of doping values in Refs. 12, O.) The most 



favorable mean-field states with 2d charge modulation 
are characterized by holes arranged in intersecting bond- 
centered stripes, i.e., the hole concentration is large on 
12 sites and zero on 4 sites, see Fig. [I]. The actual filling 
in the hole-rich regions is smaller than 1 and depends on 
microscopic parameters; due to the strong pairing cor- 
relations the 2d CDW states are good superconductors 
(with a d- wave-like pairing symmetry) . Such 4x4 plaque- 
tte states occur as large- V ground states of the Hamil- 
tonian (pi) for doping levels between approximately 13% 
and 25%7 A sample phase diagram for doping 8 — 20% is 
shown in Fig. [I] (see also Fig. 3 of Ref. |l^). For small t/J 
phase separation tendencies do not occur, therefore the 
ground state is a doped spin-Peierls state with homo- 
geneous site charge distribution. At intermediate t/J, 
frustrated phase separation leads either to stripes or to 
plaquette states; at larger t/J the kinetic energy becomes 
dominant and weakens the phase separation tendencies, 
consequently a homogeneous cJ-wave state is reached. At 
very large Coulomb repulsion, the holes arrange into an 
insulating Wigner crystal. 

The occurence of plaquette-modulated CDW states in 
favor of stripes can be understood as an interplay of ex- 
change, kinetic, and Coulomb energy as follows: The ex- 
change term prefers dimerization between spins on neigh- 
boring sites and tends to expel holes - this produces the 
overall dimer structure of the Qij, and leads to a frac- 
tion of sites being undoped. The kinetic energy is low- 
ered by possible hopping processes, i.e., by neighboring 
sites with non-zero hole density. Clearly, the exchange 
term prefers stripes, whereas the kinetic term prefers 
plaquettes, where hopping in two directions is possible. 
Now, the Coulomb energy of the 2d modulated state is 
significantly lower than that of a stripe state, because 
the charge inhomogeneity is smaller in the plaquette 
state (which is also closer to a crystalline arrangement 
of charges). From this discussion it is clear that at low 
doping, where the physics is dominated by the exchange 
term, stripe states are preferred. With increasing doping 
the kinetic energy becomes more and more important, 
which leads to 2d modulated CDW states at moderate 
values of V, before a homogeneous ci-wave superconduc- 
tor becomes the ground state. 

In all CDW states superconductivity competes with 
charge order. In the stripe states superconductivity is 
very weak due to the strong anisotropy: bulk supercon- 
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FIG. 1. Ground state phase diagram of the extended t — J 
model in the large- N limit at doping S = 20%. Thick (thin) 
lines indicate first- (second-) order transitions. All states ex- 
cept for the Wigner crystal have superconducting order. The 
stripe and spin-Peierls phases show Id charge modulation; 
the plaquette phase has charge modulations with real-space 
period 4 in both directions. The circles indicate the spatial 
distribution of hole density, the lines symbolize the strength 
of the bond variables Qij. "Full stripes" refers to states where 
the charge modulation in the large- iV limit is maximal, i.e., 
the hole density is zero in the hole-poor regions, whereas the 
"partial stripe" states have a finite hole density there. 



ductivity is established only by Cooper pair tunneling 
between the stripes. In contrast, the plaquette CDW 
states are much better superconductors, due to the full 
2d character of the charge distribution. This trend is 
consistent with the low T c in Nd-doped La2- a; Sr 2 ;Cu04 
compared to apparently charge-ordered YBa2Cu306.35 
and Bi2Sr2CaCu208+5samplesEx 



III. MODULATION IN THE LOCAL DENSITY 
OF STATES 

After having established the possible occurence of pla- 
quette CDW states in the large- TV theory for the t — J 
model, we turn to a detailed analysis of the correspond- 
ing STM signal. Hoffman et ala have introduced a STM 
technique of atomically resolved spectroscopic mapping, 
which allowed to detect LDOS modulations around vor- 
tex cores with real-space period 4, i.e., at waveyectors 
K x = (tt/2,0) and K y = (0,tt/2). Howald et alB used 
this technique to map the energy dependence, pk(w), of 
the spatial Fourier component of the LDOS at the or- 
dering wavevectors K Xj!/ . This energy dependence has 
been recently discussed within a model for the. pinning of 
SDW/CDW fluctuations by inhomegeneitiesE3, and in an 
analysis of different patterns of trapslational symmetry 
breaking in ci-wave superconductorala. n 

For a comparison with the experimental situation^ we 
restrict ourselves again to states with period-4 modula- 
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FIG. 2. Energy dependence pk(w) of the Fourier compo- 
nent of the LDOS at K = (n/2,0), obtained from diagonal- 
izing Hbcs + Hmod on a 4 x 4 unit cell. Bulk parameter 
values are t = 0.15 eV, t' = — 1/4, doping S — 17%, and a 
gap size Ao = 40 meV. The curves correspond to modula- 
tions as follows: solid - site charge density, dashed - bond 
charge density, dash-dot: pairing amplitude. The amplitude 
of the LDOS modulation is proportional to Vb in 7i mo d; here 
V<f tc = V kin = 10 meV, V r pair = 4 meV. 

tion; we will employ hopping parameters t, t' that yield 
a realistic band structure. Furthermore we have to keep 
in mind the shortcomings of the large- Y theory: the pre- 
cise location of the phase boundaries is not reliable, and 
the theory underestimates fluctuations. Therefore we will 
work with superconducting gap values close to the exper- 
imentally observed ones, and discuss both self-consistent 
mean-field solutions as well as states where translational 
symmetry breaking is imposed by hand in the Hamito- 
nian. 

To obtain initial information about the possible forms 
of Pk(w) we start by considering d-wave superconduc- 
tors with additional modulation in one of the following 
quantities: site charge density, bond charge density (ki- 
netic energy), pairing amplitude. Such an analysis 
also been independently performed by Podolsky et all 
but here we are interested in 2d modulations and further- 
more diagonalize the mean-field Hamiltonian exactly for 
the 4x4 unit cell. The Hamiltonian thus has the form 
Wbcs + ^mod, where 

Hbcs = ^2 ek4,<T C k,<r + ^2 A k(4 :T c_ kji + h.c.) (2) 

k k 

in standard notation, with Ak = Ar^cos/cz — cosfc y )/2; 
Hbcs is equivalent to the Sp(2Y) mean-field theory pre- 
sented above in the region where the large-Y ground 
state is a pure d-wave superconductor (with the cor- 



respondence Ay = Ji- 



where Aj, is the real-space 



Fourier transform of the energy gap Ak). 
The modulation is introduced via 7i mo d : 



for site CDW 



we add H s ^ d = V£ lte 



Ei /(Ri)cLcfar. fOT b0nd CDW We 



haveW^d 



^o kin E(ii) /[(Ri+RiV^c,-,, and a pair- 
ing modulation is given by H v ™ d = V£ a "J2(ij) /[(R* + 
Rj)/2](cLc}, +h.c). The function /(R) describes mod- 



ulation strength and pattern, and we will concentrate 
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FIG. 3. As in Fig. ^, but for a state with plaquette site 
charge modulation as indicated in the inset (VJf = 5 meV), 
and self-consistently determined pair fields Ay, using 
J — 70 meV. Note that LDOS modulations of similar magni- 
tude occur at wavevectors ±K X ± K„. 



on 2d bond-centered period-4 modulations with /(R) = 
[cos(7r-Rx/2 + tt/4) + cos(ttR v /2 + tt/4)]/2. 

Diagonalization of Hbcs + H mo d yields the local den- 
sity of states for each site of the unit cell, from which 
we find pvi{oS) by Fourier transformation (the real-space 
origin is chosen such that pjt{u)) is real); note that 
K = (ir/2, 0) and (0,7r/2) are equivalent with the above 
choice of /(R). Results for pk(^) are shown in Fig. || 
for the three modulation cases listed above. If we com- 
pare the curves in Fig. ||| with the STM result in Fig. 3 of 
Ref. ^ which shows a peak in the magnitude of p\t{uj) at 
subgap energies |w|/Ao 2/3, it is clear that the experi- 
ments are not well described by a site charge modulation 
alone. In contrast, our result for a modulation in the 
pairing amplitude comes closest to the curves of Ref. [|. 
We note that our results in Fig. |^ are somewhat differ- 
ent from the ones of Ref. this may be due to the 2d 
character of the modulation considered here and due to 
the approximations employed in Ref. [l5|. 

It is clear that the experimentally realized CDW state 
will have modulations in all quantities invariant under 
spin rotations and time reversal. This can - at least in 
part - be captured by a self-consistent solution of the 
mean- field equations. A natural candidate is given by 
the plaquette state found as large- TV ground state above 
(Fig. |l|). Thus, we employ the Hamiltonian Hbcs 
where the pair-field A in Hbcs (pi) is determined self- 
consistently from Aij = Jij{J 
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HJ — c ia C jph "-mod 1111 

poses a weak modulation of the site charge density as 
shown in Fig. || (the strong modulation found in the 
large- N limit will certainly be weakened by fluctuation 
corrections beyond the large- N theory). Fig. || displays a 
corresponding LDOS modulation piJ.uj). The agreement 
with the available experimental datau is not satisfying, in 
particular pk(^) does not show a large peak at energies 
below the bulk superconducting gap. 

This fact and the results in Fig. || led us to consider 
an additional effect not captured in the mean-field cal- 
culations: On general symmetry grounds, a static charge 
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FIG. 4. As in Fig. |i[ but for a model with plaquette mod- 
ulation in both site charge density and exchange interac- 
tion, and self-consistently determined pair fields Ay, using 
V sitc = 10 meV, Jo = 70 meV, V J = 3.5 meV. The dashed 
curve corresponds to a band structure with t' — — 1/3 to 
demonstrate the robustness of the result. 



modulation will lead to a real-space modulation in the ef- 
fective pairing interaction, because the CDW influences 
the local spin fluctuation spectrum. On the mean-field 
level, this can be phenomenologically accounted for by 
a modulation in the exchange interaction J. There- 
fore, we have studied self-consistent solutions of the 
mean- field theory, Hbcs + as above, but in ad- 

dition to a weak site charge modulation in H^od we 
imposed a modulation of the Jy exchange interaction, 
Jij = Jo + Vo' 7 /[(Ili + Rj)/2] for nearest neighbor sites i 
and j, which leads to corresponding modulations in both 
the pair fields Ay and the bond charge density (kinetic 
energy). Results for the LDOS Fourier component pk(^) 
are shown in Fig. |], with a rather good agreement with 
the experiments of Ref. |^. 



IV. CONCLUSIONS 

Summarizing, we have studied superconducting 
charge-ordered states of doped Mott insulators. Within 
a large- N theory we have established that ground states 
with 2d charge modulation can occur at intermediate 
doping where they are preferred over stripes. By ana- 
lyzing the energy dependence of the LDOS modulation 
as observed in STM, we have found the data of Ref. || 
to be well described by combined modulations in charge 
density as well as exchange and pairing energy, caused 
by a modulation of the pairing interaction. 
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